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We study the asymptotic behavior of the quasi-normal modes (QNMs) of w-mode pulsations of 
compact stars in the high-frequency regime. We observe that both the axial and polar w-mode 
QNMs attain similar asymptotic behaviors in spite of the fact that they are described by two 
totally different differential equation systems. We obtain robust asymptotic formulae relating w- 
mode QNMs of different polarities and different angular momenta. To explore the physical reason 
underlying such similarity, we first derive a high-frequency approximation for the polar w-mode 
oscillations to unify the descriptions for both cases. Then, we develop WKB-type analyses for them 
and quantitatively explain the observed asymptotic behaviors for polytropic stars and quark stars. 
We also point out that such asymptotic behaviors for realistic stars are strongly dependent on the 
equation of state near the stellar surface. 
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I. INTRODUCTION 

Gravitational waves (GWs) are undoubtedly a major speculation of general relativity that has not yet been observed 
directly. Decades of efforts of physicists from different disciplines have been pooled together to design and build GW 
detectors of various designs (see, e.g., [T'] and references therein). To date several Earth-based GW interferometric 
^ . detectors including GEO600, LIGO, TAMA300 and VIRGO have been operating. In spite of the fact that these 
^jT^' detectors are not yet sensitive enough to detect GWs directly, interesting upper limits have been placed on either the 
—t, GW strains or event rate for several potential astrophysical sources (see, e.g., for the results obtained from the 
latest science runs of LIGO). Optimistically speaking, GWs could be observable and analyzed directly within one or 
two decades when more advanced detectors are available. 

A promising mechanism for GW emission in astrophysics is achieved through pulsating blackholes (BHs) and 
neutron stars (NSs). When these compact stellar objects are perturbed away from their equilibrium states in various 
. astronomical processes, such as BH-BH, BH-NS and NS-NS binary mergers, and asymmetric stellar core collapse (see, 
' e.g., huge amounts of energy can be released in the form of GW, which might be detectable with our current 

(or future) GW observatories. While quantitative studies on such fierce processes demand sophisticated techniques in 
numerical relativity (see, e.g., [o'-Ts'l for some recent studies in this field), it is also popular to study the dynamics of 
the end products of these processes, usually in the form of pulsating BHs or NSs, with linearized theory, which may 
be able to throw light on analyses of numerical data gathered at the final stage of numerical simulations (see, e.g., 
for a comparison between linear and nonlinear treatments of pulsating NSs) . 
As in other linear oscillating systems, relevant physical quantities in pulsating BHs or NSs and the associated 
GW signals are assumed to have exp(— iwt) time dependence, where w is the eigenfrequency of the system under 
consideration. However, due to energy carried away by GWs, these quantities decay exponentially in time, leading to 
complex- valued frequency w = Wr + ii^i with Wi < 0. Correspondingly, these oscillation modes are termed as quasi- 
normal modes (QNMs) |17H2H . QNMs of compact stellar objects are interesting because they can, to certain extent, 
reflect the physical characteristics (e.g., mass, radius, moment of inertia, equation of state (EOS) and composition) of 
the GW emitter. In particular, a host of efforts have been paid to infer the structure of a NS, which is often masked 
by uncertainties in EOS of nuclear matter or quark matter, from its QNM spectra (see, e.g., [23 - |29| ). On the other 
hand, it has been conjectured that the asymptotic behavior of highly damped QNMs (i.e., modes with large |a;i|) of 
BHs is related to a calculation of the Bekenstein entropy in loop quantu m g ravity and to the quantum of area, thus 
spurring theorists to consider these modes with lively interest (see, e.g., [2l|, |30| and references therein). 

Motivated by the physical significance of the asymptotic behavior of BH QNMs, we study in the present paper the 
high-frequency asymptotic behavior of the (spacetime) w-mode QNMs of NSs. Although the axial and polar w-mode 
QNMs are described by two totally different differential equation systems, we find that they acquire similar asymptotic 
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behaviors and there exist robust asymptotic formulae relating these two kinds of QNMs. To theoretically understand 
these rather surprising results, we first derive a high-frequency approximation (HFA) for the polar w-mode oscillations 
to unify the descriptions for polar and axial w-mode oscillations of NSs. Subsequently, we develop WKB-type analyses 
for these QNMs and quantitatively explain the observed asymptotic behaviors for polytropic stars and quark stars. 
We also show how these asymptotic behaviors of w-mode QNMs in realistic stars could depend on the EOS near the 
stellar surface. 

The organization of this paper is summarized as follows. In Sect. II we present the observed asymptotic behavior 
of axial and polar w-mode QNMs of NSs, which is totally different from that of BHs. We also show how these two 
kinds of w-mode QNMs of NSs are related in the high-frequency regime. In Sect. Ill we introduce the HFA for polar 
w-mode oscillations and hence show that axial and polar w-mode QNMs are governed by two similar second-order 
differential systems, setting the stage for a unified treatment for the asymptotic behaviors for the two kinds of w-mode. 
In Sect. IV we establish a WKB-type analysis for relevant QNMs and derive the asymptotic behavior for w- modes of 
polytropic stars and quark stars described by the MIT bag model (see, e.g., [3l|). Through our analysis developed in 
Sect. IV, we can explain our findings summarized in Sect. II and show how such asymptotic behavior is related to the 
EOS of a star. In Sect. V we propose a symmetric and unified way to describe high-frequency axial and polar w-mode 
oscillations on equal footing. In Sect. VI we consider the case of realistic stars and discuss several key factors affecting 
our analysis. We then conclude our paper with a brief discussion and some remarks on the difference between the 
asymptotic behaviors of QNMs of BHs and NSs in Sect. VII. 

Unless otherwise stated explicitly, geometric units where G — c ~ I are adopted throughout the whole paper, and 
kilometers are used as the unit of lengths. Besides, we will only consider non-rotating NSs (with ideal fluids) and BHs 
in the present paper. As usual, the geometry of a static gravitational source with spherical symmetry (NS or BH) is 
described by the line element: 

ds^ = g^pdx^dx^ = -e^dt^ + e^dr^ + r^{dd^ + sin^ dd(f>^). (1.1) 

For NSs the metric coefficients e'^^'"' and e"^^''-' can be obtained from the equilibrium configuration governed by the 
Tolman-Oppenheimer-Volkov (TOV) equations [32l. Issj. 



II. ASYMPTOTIC BEHAVIOR OF W-MODE QNMS 



The oscillation modes of both NSs and BHs can be classified into polar and axial classes according to the property 
of the metric perturbations under inversion. For typical stellar models, the polar class includes (fundamental) f-mode, 
(gravity) g-mode, (pressure) p-mode and (spacetime) w-mode (see, e.g., for a review). Unlike the f-, g- and 
p-modes which have their Newtonian origins 34], the w-modes mainly represent oscillations in spacetime and are 
specific to theory of general relativity. In contrast to the rich spectra of polar oscillations, axial oscillations of NSs do 
not introduce non-trivial motion of matter and all are classified as the w-mode. 

Axial QNMs of NSs are governed by a second-order differential equation, namely the Regge- Wheeler equation for 
NSs [3i,[3^: 



drl 



V'a(»^*) = 0, 



where the tortoise radial coordinate r, is related to the circumferential radius r through the relationship: 



(2.1) 



(2.2) 



and the NS Regge- Wheeler potential is given by: 



^ [l{l + l)r + iTTr^{p-p) - 6m] 



(r - 2M) 



[l{l + l)r - 6M] , 



< R: 



r> R. 



(2.3) 



Here M, i?, p{r) and p{r) are the total mass, the radius, the density and the pressure of the star, / is the angular 
momentum index, and 



m(r) = An p{r')r' dr' . 



(2.4) 
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Note that for r > R the NS Regge- Wheeler potential is identical to that of the BH case. 

In contrast to the simplicity of axial oscillations, polar oscillations of NSs are more complicated as oscillations 
in spacetime and fluid are in general coupled together. A set of coupled fifth-order equations gover ning non-radial 
pulsations of compact stars were first derived in 1967 by Thorne and Campolattaro |37| . Based on [37[, Lindblom 
and Detweiler reduced the system into a fourth-order one fssl Isol. Subsequently, there have been several alternative 
formulations for polar oscillations of compact stars (see, e.g.. [36l l45 - l45l |). and most of them (if not all) are fourth-order 
systems. 

As introduced above, axial and polar oscillations are governed by two completely different formalisms. Therefore, 
it seems plausible that axial and polar w-mode QNMs would acquire different asymptotic behaviors. We have carried 
out numerical calculations to locate these QNMs for various stars, including polytropic stars (with EOS p oc p^), 
quark stars (with EOS p = (p — 4i3)/3) and NSs constructed with several popular EOSs of nuclear matter, to examine 
relevant high-frequency behavior. The results obtained are in fact a bit surprising and encouraging as well. As shown 
in Fig. [l] where axial and polar w-mode QNMs of quark stars and polytropic stars are shown and analyzed, for a 
specific star these two kinds of w-mode QNMs lie on a common smooth curve star in the complex w-plane, which can 
be approximated by the following equations: 

(2.5) 

oj-A = i In [rfc(i?*)-^+2] - (y + 1) H^rR*), (2.6) 

where j — 1,2,3,... is the mode number labeling QNMs in order of increasing frequency, 9i, and > are real 
parameters dependent on the star, a = 0(1) for axial (polar) w-mode, TV = 1/(7 — 1) is the so-called polytropic index 
{M = for quark stars hereafter), and i?* is the tortoise radius of the star. 

The situation for realistic stars is less uniform and it is difficult to extract a generic asymptotic behavior for Wj (see 
Fig. 121) ■ However, uj^ still obeys an approximate equation similar to (12. 5p . namely 

(2.7) 

Here the parameter i?* is usually a bit less than and could have a mild dependence on uj. Despite the lack of 
an explicit formula for cjj, there is a universal relationship between the complex QNM frequencies of axial and polar 
w- modes: 

.['^ = {.^^^+.^^1)/^ (2.8) 
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and vice versa, where the superscripts P and A denote polar and axial modes, respectively. This equation holds 
asymptotically for large j and works nicely irrespective of EOS of the star (including various nuclear matter EOSs, 
polytropic EOS and quark matter EOS). The main objective of the present paper is to establish a proper theory to 
explain these observed asymptotic behaviors, which will be the focus of the following sections. In fact, we will show 
how the values of 9b and rt for polytropic stars and quark stars can be obtained analytically. 



III. WAVE EQUATION FOR POLAR W-MODE 



The similarity between the asymptotic behaviors of the two kinds of w-mode is a bit surprising since they are in 
general governed by two entirely different sets of equations. The findings discussed in the last section strongly hint 
that axial and polar w-modes could actually be placed on an equal footing at least in the large- frequency limit. In this 
section we shall develop a HFA for polar w-modes under the inverse-Cowling approximation (ICA) [4^ \^ to show 
that they can be described by a single second-order wave equation analogous to (|2.ip . To this end, we first introduce 
an alternative formalism of non-radial polar oscillations of compact stars suggested by Allen, Andersson, Kokkotas 
and Schutz [iJl, referred to as AAKS formalism in the following discussion. 
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A. AAKS Formalism 

In AAKS formalism polar oscillations of compact stars are described by two spacetime metric variables F, S and 
one fluid perturbation variable H: 



rF = hgg, 

rS = htt Trha 



H = 



Sp 

P + p' 



(3.1) 
(3.2) 

(3.3) 



where haeia, P = t,r, 9, (jj) denotes the first-order perturbation in the metric, and Sp is the Eulerian variation in 
pressure 44']. These three physical quantities satisfy three second-order evolution equations plus one second-order 
time-independent (or equivalently frequency-independent) Hamiltonian constraint equation and one can use any three 
of these four equations to solve for F, S and H. In this paper we shall use the two wave equations: 



(PS 26" 



S+-rY + — [27rr3(p + 3p) + m - {n + l)r] S - 



4e' 



2u 



drl 
dr'i 



3\2 



(m -|- Airpr^) 
r — 2m 



Anpr^ — 3™ 



uj^F+^^ + — [2T:r^{3p + p) + m- {n + l)r] F = -2 [AT:r^{p + p) - e-^] S + 8Tr{p + p^e" ( 1 



(3.4) 



(3.5) 



plus the Hamiltonian constraint: 

PF 



,(i/+A)/2 



dF 



= 5 {m + AT:r^p)- \- —[Uirr^ p - m - 2{n + l)r]F 

art ar^ r-^ 



.2 

dS_ 



8TTr^{p + p) -{n + 'S) + 



4m 



S+^e'^{p + p)H 



(3.6) 



to develop a second-order wave equation for polar w-mode. In these equations n — {I — + 2)/2, and Cs is the 
sound speed in the stellar fluid. 



B. Inverse-Cowling approximation 



It is commonly accepted that polar w-mode oscillations are oscillations in spacetime and fluid elements are to 
certain extent spectators. Hence, it is reasonable to expect that polar w-mode QNMs could be reproduced under the 
so-called ICA where certain fluid perturbations (e.g. fluid displacements or pressure variation) are neglected (icl. |47|. 
In fact, it has been shown that ICA assuming zero Lagrangian pressure variation can produce quantitatively correct 
result 47]. Here we propose another version of ICA, namely zero Eulerian pressure variation (i.e., = 0), to derive 
a valid wave equation for polar w-mode QNMs. 

From p.5|) and p.6p . we have: 



Snre^Sp = cj^F + (m + iTrpr^)— + —(3m -|- 'inpr^)F + re-^^+^^/^-i + {n + 1)S. 

dr^ dr^ 

Under our ICA scheme, Sp — 0, and consequently Eq. p.7p can be rewritten in the following way: 



F 



ei^+>^)/^(m + 4npr^) dF e" {im + Airpr^) re-^^+^^Z^ dS n + 1 

F ; — S. 



This equation, together with (13.41) . can lead to a single wave equation in S. 



dr^ 



(3.7) 



(3.8) 



C. High-frequency approximation 



We note that the frequencies of w-mode are at least a few times higher than that of /-mode, which is proportional 
to (M/i?^)^/^ ~ P"'^^^: with p being the average density of the star. As a result, terms of the order p^/^/w are small for 
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w-mode QNMs and can be neglected for high-order modes. Based on this observation and (I3.8p . we have developed 
a systematic power expansion in l/w for F. Hence, the contribution of F to (j3.4p can be included perturbatively, 
leading to an accurate scheme to solve for S. The details of such perturbation scheme are provided in Appendix A 
for reference. However, since we are interested in the high-frequency asymptotic behavior of polar w-mode, we adopt 
the lowest order HFA where the contribution of F to (|3.4p is neglected completely to obtain a single wave equation 
for S: 



5(rO = 



(3.9) 



with 



[l{l + l)r-ATTr^{p + ip)-2m\, r < R; 
+ l)r - 2M] , r>R. 



(3.10) 



This polar w-mode wave equation is identical to the wave equation for the axial mode, (|2.ip . except for the functional 
form of the potential term. We will see in the following discussion that such similarity between the governing equations 
of these two kinds of w-mode wave is in fact the origin of the peculiar asymptotic behavior observed in Sect. H. 

To gauge the accuracy of HFA and examine the validity of p.9|) . we located the first ten polar w-mode QNMs 
with (|3.9I) and compared the numerical results with the exact values in Tables U and HI] for a polytropic star (EOS: 
p — 100/9^) with a compactness 0.297, and a realistic star (EOS: APR2 [4^) with a compactness 0.196, respectively. 
The numerical results of u are given in units of M^^ . From the numerical results, it is obvious that the error of HFA 
decreases with increasing oj. Therefore, the HFA can indeed provide us accurate results for high-order modes. 



IV. ASYMPTOTIC BEHAVIOR OF W-MODES FOR MODEL STARS 

In the above discussion we have shown that both axial and polar w-mode QNMs of compact stars are governed by 
a Klein-Gordon type wave equation: 



^(r*) = 0, 



(4.1) 



where V — (Vp) and ip — ijjs. [S] for axial (polar) case. Despite the fact that the form of (14.11) is similar to the 
Regge- Wheeler equation and the Zerilli equation describing GW propagation around BHs [H, |49{, the potential V has 
non-trivial dependence on density and pressure. To see how such dependence could affect the asymptotic behavior 
of relevant QNMs, in this section we focus our attention on stars with simple analytic EOS, namely polytropic stars 
and quark stars described by the simplest MIT bag model with zero strange quark mass. 



A. Locating QNMs 

One usual way to locate QNMs is as follows. For any specific oj, we first find a function /(aj,r,) satisfying the 
regularity condition / ~ r^"*"^ as r* — ?■ and a function g(w,r*) satisfying the outgoing condition g ^ exp(iwr*) as 
r* oo. QNMs are then located when the Wronskian of f{uj,r^,) and g{uj,rt) vanishes, i.e., 

W{f,g;uj) = f{uj,n)-^g{uj,r,) - g{uj,r,)-^f{uj,r,) = 0, (4.2) 

and the corresponding eigenfunction is proportional to f{u!,r^) (or equivalently g(w,r*)). 

In the high-frequency domain the asymptotic expression of f{uj,r^,), g{uj,r^,) and W{uj) can be derived under the 
standard WKB approximation, and hence an approximate expression of QNM frequencies can be derived as well 
[sol [sTj . As in other applications of the WKB approximation, due attention has to be placed on the singular points of 
the potential, e.g., points at which the potential (or the derivatives of the potential) are discontinuous. For polytropic 
stars and quark stars, V{rf,) is well behaved except at the stellar surface. We shall show in the following analysis how 
the reflection of GW at the stellar surface could dominate the asymptotic behavior of w-mode QNMs. 
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B. WKB analysis 



Adopting the methods discussed in [52, [SJl , we have developed a systematic WKB analysis for both / and g (see 
[5^ for details). In the high-frequency limit, the dominant parts in / and g are given respectively by 



/(r,,tj 00) ^ uj7\ji{ujr^), 



< < i?* 



and 



g{r^,,uj 00) = 



exp [iw(r, — i?,)] + 7?.exp [— iw(r, — i?*)] , l/uj ^ r» < i?*; 



(1 + TZ) exp [icj(r, — i?,)] 



-R* < r, < 00. 



(4.3) 



(4.4) 



Here 7?. is the standard reflection coefficient as defined in [52|. In particular, if the potential V{r^,) has a discontinuity 
in d^V/dr^ {N = 0, 1, 2, . . .), 7^ = aw-(^+2), with 



^+2 



lim 



dr^ 



\r,=R,+e 



(4.5) 



Detailed derivation of the above result can be found in [5^ . 

The asymptotic eigenvalue equation for w-mode QNMs is then obtained by setting the Wronskian of ^(r^, w 00) 
and /(r* , lu 00) to zero: 



= exp(-i2wi?,), 



(4.6) 



where h = ri,exp(i6i,) = (— l)'+^a with > and 9), being two real constants. This equation can be solved readily, 
yielding 



27rj — 9b — i 



(4.7) 



where j 3> 1 is a positive integer and the approximation |aj| ~ \uj^\ has been used. 

It is obvious that Eq. (14. 7p indeed leads to the asymptotic behavior of w-mode summarized in Sect. II as long as N 
is equal to the polytropic index TV. In the following discussion we shall show that this is the case for A/" = 0, 1, 2, . . .. 
More importantly, Eq. (|4.7p can be generalized to include cases with non-integral values of TV. Before diving into 
detailed calculation of 6, several useful features of (|4.7p are worthy of remark. Firstly, uj-c of the j-th and (j -f l)-th 
eigenfrequencies are spaced by tt/R^. Secondly, depends linearly on the logarithm of The above features both 
hold for axial and polar high-frequency w-modes. In addition, since the leading discontinuity in the potential or 
V-p is usually caused by the ±ATie^ p term, the expression of h has a minus sign difference between the axial and the 
polar cases. The phases of b for the axial and polar cases thus differ by tt, leading to a shift of 7r/(2_R») in oji-. All 
these features are consistent with our numerical results presented in Sect. II. 



C. Application to model stars 



To make our discussion concrete and physical, we apply our analysis to quark stars and polytropic stars, which are 
described by EOSs with analytic forms. It is readily shown that the potential and the derivatives of the potential for 
such stars are continuous except at the stellar surface. In the following we will show how the values of N and b can 
be expressed in terms of the physical parameters of the star. 

For a quark star with the MIT Bag Model EOS (p = 3p + 4i?, with B being the bag constant), the star surface 
is defined at the point where the pressure vanishes. However, the density there is equal to 4i? and does not vanish. 
Therefore, the potential has a direct discontinuity at the surface of the star, i.e., the = limit of ()4.5p . and 
straightforward calculation shows that 

b = {-ly+'^ABiT (^1 - ^ 

We now turn our attentions to polytropic stars with EOS given explicitly hy p — Kp^^^l^ ^ where k is a constant. 
Again, it is easy to see from TOV equations that the variables m, p, -p and v of polytropic stars are all analytic 
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functions inside the star and hence our theory apphes. Near the surface of the star, following directly from the TOV 
equations, the density and the pressure behave respectively as (38| : 



P 
P 



where 



M 



A/" 



(4.9) 
(4.10) 



(4.11) 
(4.12) 



We first consider a simple case where A/" is equal to a positive integer N. Then the lowest order discontinuity appears 
in the iV-th derivative of the potential at the surface and the coefhcient b can be found explicitly: 



6=(-l)'+-vrr(A/-+l) (l-^) 



-iM 



where r(2;) is the gamma function. It is interesting to note that Eq. (j4.13l) in fact includes 
A/'= 0. 

It follows directly from the asymptotic analysis developed above that for j 3> 1 



LO-^R^, = 



(y + l)ln(c^,i?,) 



■ In [rb{R, 



\^f+2-\ 



with 



(TTT 



-I 



2 



(4.13) 

as the special case 

(4.14) 
(4.15) 

(4.16) 



It is obvious that the forms of (|4.14p and (|4.15p agree with the numerical results summarized in Sect. II. Most 
interestingly, Eqs. (|4.13p . (|4.14p and (I4.15P can be readily generalized to non- integral values of TV by method of 
analytic continuation. Thus, we have established a generic and unified asymptotic behavior for polytropic stars and 
quark stars as well. 



D. Numerical results 



To further examine the validity and accuracy of our method, we have calculated the exact eigenfrequencies for 
w-mode QNMs of several model stars (including one quark star and six polytropic stars whose physical parameters 
are tabulated in Table IIIII for reference) and showed the results in Fig. [T] It is clearly shown in these figures that 
the asymptotic behavior of w-mode QNMs agrees nicely with the theoretical results derived above. For example, the 
spacing in uj^ between two successive QNMs of the same polarity is close to ir/R^,. There is one axial (polar) QNM at 
the midpoint of the line joining two neighboring polar (axial) QNMs (i.e. Eq. p.SI) ). Besides, Wi depends linearly on 
In j. These features are just the direct consequence of (|4.14p and (|4.15p . 

On the other hand, the coefficient b for each case can be calculated from (|4.13p with the physical parameters 
tabulated in Table Hill Substituting the so-obtained value of b (see Table Hill) into (|4.14p and (|4.15p . we can determine 
both Wr and Wj for j ^ 1 without using any unknown free paramter. To further verify our theory, we have fitted exact 
values of Wr and Wi with the following functions: 

uJrR* = cij + Co, (4.17) 
Wii?* = dilTi{ujrR*) + do, (4.18) 

with ci, Co, di and do being adjustable parameters to yield the best fit. The theoretical values of these parameters 
(see (|4.13l) . (|4.14p and (14.151) ) and the values obtained from the best fit to the w- modes (from the 11-th to the 20-th 
modes) are listed in Tables ITVl (values of ci, cq) andlVl (values of di, do). It is obvious that our analysis works nicely 
in all cases. The numerical values of both ci and di are close to their theoretical ones. For the case of co and do, the 
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errors are larger, especially for the former. This is expected as sraall errors in the slope would result in much larger 
errors in the intercept in the process of extrapolation. Therefore, as a remedial action, we have first set ci and di to 
their theoretical values and substituted the eigenfrequency of the 20-th w-mode into (|4.17|) and ()4.18|) to compute cq 
and do directly (the numeric results enclosed in brackets in Tables IIVI and |V| . It is clear that this modified method 
can improve the accuracy to certain extent. 



V. UNIFIED APPROACH TO W-MODES 



From the asymptotic analysis on w-modes developed above we see that the asymptotic behaviors of w-modes 
are mainly determined by the most singular part in the potential, which is responsible for refiection of GWs at 
high frequencies. Apart from the trivial centrifugal barrier that has been duly considered by adopting the solution 
/ oc j{Lor^), there are three other terms in the potential which are proportional to p, p and m, respectively. As usual, 
the effect of p is much less than that of p. Moreover, m is an integral over the p term (see (12.41) ) and should be less 
singular than p. Hence, it is reasonable to argue that both the contributions of p and m to the potential are small 
at high frequencies. Once the terms proportional to p and m are neglected in the potential, 14 and V-p are almost 
identical except for a sign difference in the term proportional to p. More importantly, such argument should work 
irrespective of the EOS of a star. We therefore propose to use the following unified approach to w-modes: 



■0(w,r*) = 0, 



(5.1) 



with 



+ l)r + {-ly ^nr^ p] 



-l{l + l){r-2M), 



r <R- 



r> R. 



(5.2) 



Thus, the two kinds of w-modes are cast in a symmetrical description. 

We have tested the validity of our proposal by comparing the numerical results of the w-mode QNM eigenfrequencies 
obtained from the above approximate wave equation with the exact values in Tables I VII (axial case) and IVIII (polar 
case). The star under consideration is a realistic NS with EOS A [55|, which has a compactness of 0.27 and a central 
density pc = 2.227 x 10~^. The w-modes listed in the tables are the 3-rd to 12-th modes. One can see that the unified 
approach works accurately for both kinds of w-mode. The percentage error in uj-c and Wj is less than 1% in most cases 
and decreases with increasing order. Thus, the validity of the approximate form of the potential in (|5.2I) is justified. 
As a result, w-modes of both polarities can be cast into a unified description. 



VI. GENERALIZATION TO REALISTIC NEUTRON STARS 



A. Numerical result 



It is not easy to generalize the method established above to analyze the w-modes of NSs characterized by realistic 
EOSs of nuclear matter for the following reasons. First, such EOSs are usually given in the form of discrete data 
points and methods of interpolation are needed in order to generate the function p{p) (see, e.g., [1^). It is likely that 
artificial discontinuities in the potential (or its derivatives) may develop when p{p) so generated is used as the input 
to solve the TOV equations describing the static configuration of a compact star. These artificial discontinuities could 
then lead to extra reflections of GW and obscure numerical results of high-order QNMs and the related theoretical 
study as well. Second, the polytropic index of nuclear matter is in general density-dependent. However, as shown in 
the previous section, the polytropic index is a crucial factor in the asymptotic behavior of w-mode QNM. Except for 
GWs with very high frequencies, reflections occurring at interfaces other than the stellar surface are characterized by 
different polytropic indices, leading to more uncertainties in the asymptotic behavior of realistic NSs. 

Instead of deriving the asymptotic behaviors of w-modes for realistic compact stars from first principles, we have 
gathered some empirical rules for high-frequency w-mode QNMs based on numerical studies carried out for numerous 
stellar models, including models A 55 1, C [l^l, three models (AU, UU and UT) proposed in [sl], and models APRl 
and APR2 proposed in '48']. In particular, we have constructed NSs based on the two analytic EOS proposed by 
Haensel and Potekhin j5^], which can well approximate EPS EOS [gO] and SLy EOS ^] of realistic neutron-star 
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matter. More importantly, these analytic EOS are free from artificial discontinuities due to interpolation of discrete 
data points and can facilitate theoretical analysis. As a representative example to illustrate our discussion, we show 
in Fig.[2]the w-mode QNMs of two NSs obeying the analytic SLy EOS with compactness C = 0.306, 0.228, from which 
one can see the following features. First, despite the fact that Wr no longer follows (|4.14[) . it obeys a modified version 
of (|4.14p . namely (j2.7p . where i?* is replaced with R„. Strictly speaki ng, i? ^ also has a mild frequency dependence, but 
such dependence is not obvious in the frequency range shown in Fig. 2(a) Second, uj{R^ is no longer a linear function 
of ln(ajiJ?*). The graph jwii?, | versus ln(a;r-R*) shown in Fig. 2(b) displays a frequency-dependent slope. The slope of 
the graph increases from about 1.3 (1.3) to about 1.4 (1.9) for the case C — 0.306 (C — 0.228) when j goes from 5 to 
20. Generally speaking, the slope demonstrates a greater dependence on frequency for NSs with smaller compactness. 
Third, although for realistic stars both lj^ and cjj do not follow (|4.14l) and (|4.15p . the asymptotic behavior Eq. (|2.8p 
still holds. Hence, axial and polar w-modes are located on the same smooth curve and evenly alternate with each 
other in the high-frequency regime. In the following we try to explain the above findings in a semi-quantitative way 
based on the asymptotic analysis developed for polytropic stars. 



B. Physical picture 

First of all, it is instructive to associate such asymptotic behavior of w-mode QNMs with the physical structure 
and the EOS of realistic NSs. As is well known, a typical realistic NS consists of a core with supernuclear densities 
(i.e. p > 2.8 X lO^^gmcm^'^), and a thin crust at subnuclear densities. The crust can be further divided into an inner 
part and an outer part by the neutron drip density p « 4 x 10^^ gmcm"'^. Due to the effect of neutron drip, stellar 
matter softens significantly and the adiabatic index F, a typical measure of stiffness defined by 

r^(i + ^)^^, (6.1) 

V Pj Pdp 

decreases abruptly from about 1.4 to about 0.4 near the neutron drip density. On the other hand, at the crust-core 
interface nuclei disappear and nuclear matter there consequently becomes particularly stiff owing to the repulsive part 
of short range nuclear force. Correspondingly, F experiences a sharp rise from about 1.4 to about 3 around the density 
p ^ 2 X lO^^gmcm^'^. These softening and stiffening of nuclear matter give rise to huge and rapid changes in the 
density profile in the inner crust, which typically has a thickness of only about a few percent of the stellar radius. For 
GWs with wavelengths greater than the thickness of the inner crust, these drastic changes in the density effectively 
lead to discontinuities in the scattering potential V and its derivatives as well. Hence, strong reflections of GWs are 
expected to occur at the core-crust interface and such a mechanism indeed attributes to formation of QNMs. 



C. Semi-quantitative analysis 

As mentioned above, for a realistic NS, GWs are reflected not only at the stellar surface but also from other abrupt 
changes in p inside the star. In general, the reflection coefficient TZ in ()4.4p can be obtained from the leading term in 
Eq. (10) of Ref. [H, which yields: 

n = ^"P(f!"^*) r dr. ^ exp(2^.r.) , (6.2) 

where in V{r) the trivial centrifugal barrier has been excluded. This equation is generally valid and reduces to (|4.5|) 
if V (or its 7V-th order derivative) has a step discontinuity at r. = R.. It can be clearly seen from (|5.2p that the 
potential V is dominated by the density term p. Hence, the reflection coefficient TZ is proportional to the Fourier 
transform of dp/ dr.. 

In Fig. [3] we show dp/ dr. versus r. for the two NSs constructed with the analytic SLy EOS [5§|. It is interesting 
to see that there is a sharp peak in \dp/dr.\ at r. = Rm near the surface of the star. Generally speaking, Rm/R* 
increases slowly while the width of the peak decreases with increasing compactness. This peak obviously reflects 
the abrupt change in the stiffness of nuclear matter around the core-crust interface. In fact, it can be argued that 
Rjn = i?2, where at r. = i?2 the adiabatic index F equals 2 (see Appendix B and the data shown in Table IVlTl|) . 

For QNMs with wavelengths much longer than the width of the peak, dp/dr. can be well approximated by a S- 
function. Following directly from standard theory of asymptotic expansion of Fourier integrals (see, e.g., [62)]), the 
leading behavior of the reflection coefficient is 

TZ ^ exp [2iuj{R^ - R*)] , (6.3) 
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where ao(w) approaches a constant at low frequencies and decreases rapidly if l/|a;| is much smaller than the width 
of the peak in \dp/dr^\. Apart from this leading term, due to the asymmetry in \dp/dr^:\ around r» = i?™, there are 
other sub-dominant terms in the the reflection coefficient, leading to the following asymptotic of TZ: 

7e « f ^ + ^ + ^ + . . .) exp [2iLo{R„, - i?,)] • (6.4) 

Physically speaking, the coefficients ao(w), ai(w) and 02(0;) measure the discontinuities in p, dp/dr^, and d'^p/dr^ 
around the point r* = Rm, respectively. As in the case of ao(a;), ai(w) and 02(0;) decreases rapidly if l/\uj\ is much 
smaller than the typical length scale associated with the corresponding discontinuity. 



D. Asymptotic behavior 

Recalling the asymptotic eigenvalue equation (j4.6|) for polytropic stars and taking (j6.4p into consideration, we can 
show that for realistic stars the w-mode QNMs are asymptotically given by the solution of the equation: 

n{uj) w (-1)'+^ exp(-i2a;i?„). (6.5) 

Since the most singular part of the potential is dominated by the p term, from (|4.5|) and (|5.2p we see that the 
generalized reflection coefficients for axial and polar w-modes should have the same frequency dependence but with 
a sign difference. Therefore, the two eigenvalue equations for axial and polar w-modes can be unified into one single 
equation: 

TZ^iiu) « (-1)'+!+- exp(-^2cJ,i^,„) exp{2iUiR„,) , (6.6) 

with Tin being the reflection coefficient of the axial w-mode, and cr = 0(— 1) for the axial (polar) case. In the 
asymptotic region, T^a is further assumed to have a constant phase, i.e., TZs.{uj) « exp(i6'a)|7?.a(w)|, Eq. (|6.6p then 
leads to: 

uj,R^ = TTj , (6.7) 

uJiRm = iln(|7^a(c^)|). (6.8) 
Upon equating ^a ^ (' + l)?!" with 9b, Eq. (j6.7p agrees nicely with (j2.7p . which is an empirical relation obtained 



from numerical data such as those shown in Fig. 2(a) provided that i?, = i?„j. Our argument readily provides the 
numerical value of i?* in (12.71) a proper physical interpretation, i?* is actually the tortoise coordinate of the point 
at which \dp/dr^\ attains a maximum. In addition, i?, is also close to the tortoise coordinate of the point where the 
adiabatic index is 2 and nuclear matter becomes stiff due to the disappearance of individual nuclei. This point is 
conflrmed by the numerical data shown in Table IVIIII 

As suggested in (|6.8p . Wi depends crucially on the magnitude of the reflection coefficient Tia.{uj). For the cases of 
polytropic stars and quark stars, 7^a(w) oc 0;^''^+^', leading directly to (|4.15p . and the plot of versus ln(wr-R*) 

is a straight line with a slope rus = [N + 2)/2. In the cases of realistic NS, the polytropic index M as well as the 
adiabatic index V depend on density p. It is therefore reasonable to expect that the plot of ui versus In Wr is no longer 
a straight line. This is in agreement with the result shown in Fig. |2(b)[ As a matter of fact, the slope of the curve, 
ms(ajr) = [N'mi'^r) + 2]/2, increases gradually with Wr and the trend is more obvious for the star with C — 0.228. In 
Fig. Uthe value of Nm for the is plotted against ln(wi.i?*) two SLy stars. As shown in the following discussion. Mm 
so obtained could indeed measure the variation of the polytropic index M near the stellar surface. 

Since the reflection of GWs mainly takes place around the stellar surface, we define the average value of F over a 
surface layer with a thickness measured by the wavelength A = 271 /uj^, namely: 



1 rR. 

r=-/ r(r,)dr., (6.9) 

^ -'R.-A 



where A = rjX with rj = 0.61 (0.7) for the star with C — 0.306 (0.228) in our calculations. In Fig. |3] we compare the 
values of Afm with 77 defined by 



(6.10) 
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It is interesting to note that can indeed approximate Mm nicely. Thus, the hnk between the asymptotic behavior 
of w-mode QNMs and the EOS of the nuclear matter distributed within around one wavelength about the stellar 
surface is clearly displayed. 

On the other hand, the asymptotic expansion of the reflection coefficient in (|6.4p could also cast light on the 
variation of the slope m(a;r). For GWs with wavelengths longer than the thickness of the crust, they get reflected 
from the narrow peak in \dp/dr^\ around the core-crust interface. Such reflection is measured by the term ao/w'^ in 
the expansion of TZ. However, as shown in (|6.4|) . there are other contributions ai/w^ and 02/10^ due to reflections 
resulting from discontinuities of higher orders. Each of these terms alone leads to an effective slope 1, 1.5 and 2, 
respectively. On the other hand, for GWs with wavelengths less than the width of the peak, the value of ao decreases 



markedly. Thus, we expect that the frequency-dependent slope of the graph versus ln(wr-R*) in Fig. 2(b) 

increases with frequency gradually from 1 to 1.5, and then to 2. Such increase in the slope is expected to be more 
pronounced if the peak in dp/dr^ is broader. Our conjecture indeed agrees with the numerical data qualitatively. 
Furthermore, we deflne an effective (frequency-dependent) polytropic index Me by: 

M^-'-^-2. (6.11) 

O In CJr 

For polytropic stars Me reduces to the standard polytropic index TV. As shown in Fig. 21 the trend of Ne (minor 
oscillations in Ne have been suppressed) agrees qualitatively with those of Mm and M. This once again demonstrates 
the close relationship between the high-frequency behavior of w-mode QNMs and the mass distribution of a NS. 



E. Methods relating QNMs 

In spite of the complexity in the situation of realistic NSs, the relationship between the axial and polar w-modes for 
these stars is identical to that for polytropic stars. Due to the sign difference in the reflection coefficients of axial and 
polar GWs, lu-^ of axial and polar w-mode frequencies, with the same mode index j, differ approximately from each 
other by 7r/(2i?,). Hence, the real part of a polar w-mode eigenfrequency lies in the middle of two axial w-modes and 
vice versa. The difference of Wr leads to the difference in through (16. 8p . Therefore, the imaginary part of a polar 
w-mode eigenfrequency should also lies approximately in the middle of two axial w-modes. 

Based on the relationship between the two types of w-mode, two methods are proposed here to derive approxi- 
mate polar w-mode eigenfrequencies from the corresponding axial w-mode eigenfrequencies. The first method is to 
approximate u> of the j-th polar w-mode QNM with the mean value of the eigenfrequencies of its two neighboring 
axial w-mode QNMs of the same I, i.e., wj'^' = ('^j"*' +wj+i)/2. Besides, from (|6.7p it is obvious that the asymptotic 
behaviors of polar w-modes with I = Iq and axial w-modes with I ^ Iq + 1 are the same. Hence, it is possible to 
approximate the frequency of a polar w-mode with I = Iq with that of an axial w-mode with the same mode order 
but I = Iq -\- 1. As numerical calculations of axial w-mode eigenfrequencies are usually much easier and faster than 
polar w-mode calculations, the two methods mentioned above are deemed useful to provide quick ways to locate 
approximately polar w-mode QNMs. As an example to illustrate our method, we consider the w-mode QNMs of a 
model C NS (with compactness C — 0.186 and central density pc = 8.752 x 10^'') and evaluate approximate polar 
w-mode frequencies (from the 5-th to 14-th modes) for the case Iq — 2 using these two methods. The approximate 
results so obtained are listed and compared with the exact values in Table IIXI The accuracies of the data shown in 
the table are impressive especially for high-frequency modes. 



VII. CONCLUSION AND DISCUSSION 



We study the high-frequency asymptotic behavior of w-model QNMs of compact stars, including model polytropic 
stars, quark stars described by the simplest MIT bag model and NSs constructed with various realistic EOSs of 
nuclear matter. By examining relevant numerical data, we find that there exist some model-independent generic 
behaviors in these QNMs. Most interestingly, notwithstanding the difference between the formalisms of axial and 
polar w-modes, the asymptotic behaviors of QNMs of these two kinds of w-mode are quite similar. In fact, as shown 
in (12. 8|) . the spectrum of polar w-modes is derivable from its axial counterpart and vice versa. Such coincidence has 
motivated us to study the physical link between these two kinds of w-mode QNMs. Based on the AAKS formalism 
for polar oscillations of compact stars [4J|, we propose another version of ICA (i.e. H = 0) and introduce HFA (i.e. 
F = 0) to derive a Klein-Gordon type wave equation p.9p for high-frequency polar w-modes. Moreover, noting the 
physical significance of the density term in the scattering potential, we show that the potentials in the axial wave 
equation (j2.1[) and the HFA polar wave equation (j3.9|) can be unified and rewritten into similar forms as in (15.21) in 
the high-frequency regime. 
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Having established a unified approach to both axial and polar w-modes, we formulate a WKB analysis to evaluate 
the eigenfrequency for relevant QNMs. We note that the crucial factor affecting the analysis is the discontinuities 
in the potential (or the derivatives of the potential) . We then apply our theory to the cases of polytropic stars and 
quark stars which are described by smooth EOS. Simple yet accurate formulae (|4.14p and (I4.15p . respectively for Wj- 
and Wi, are obtained. Despite the fact that similar development for realistic NSs is less trivial and often plagued 
by the artifacts in EOS, we can still provide a qualitative explanation for the asymptotic behavior in such cases. In 
particular, we point out the relationships among the asymptotic behavior, the polytropic index of nuclear matter and 
the Fourier transform of dpjdr^,. 

It is worthy to note that the asymptotic behavior of QNMs could easily be masked by numerical errors in the 
density profile. Such errors could lead to artificial reflections of GW and hence unexpected features in the distribution 
of high-frequency QNMs. In order to accurately locate high-frequency QNMs, due care has to be taken in solving 
the TOV equations for the equilibrium configuration of a star (see [5J] for a detailed study for the effect of a tiny 
truncation error in the density distribution near the stellar surface on the asymptotic behavior of QNMs of NSs). 

Lastly, we would like to compare the asymptotic behaviors of QNMs of NSs and BHs. Unlike oscillations of NSs, 
which can be classified into fluid and spacetime modes, oscillations of BHs are all spacetime modes. Besides, axial 
and polar mode oscillations of BH are governed by two different equations, namely the Regge- Wheeler equation and 
the Zerilli equation [H, H^. However, it is well known that the spectra of these two kinds of QNM are identical due 
to the symmetry in the scattering potentials [63l - [65j . In fact, the potentials in the Regge- Wheeler equation and the 
Zerilli equation are super-symmetric partners of each other 65]. However, such symmetry between the axial and polar 
w-modes of NSs no longer exists. Therefore, it is interesting to note that our finding in the present paper partially 
restores the symmetry between these two kinds of QNMs in the sense that their spectra are closely related through 
(12. 8p in the high-frequency regime. On the other hand, for QNMs of BHs the frequency of the j-th QNM (j ^ 1) is 
asymptotically given by (see, e.g. [lllls^l and references therein): 

SttA/wj w ln3 - + 1/2) + 0(1/V7) ■ (7.1) 

Therefore, is bounded and approaches a constant, while Wi is proportional to j + 1/2 and can increase indefinitely. 
This is in sharp contrast to ()4.14|) and ()4.15|) . the formulae for NS w-mode QNMs. Such distinction between QNMs 
of BHs and NSs is attributable to the analytic property of relevant potentials in the wave equation. For the BH case 
the potential is smooth everywhere for oo > r > 2Af and gravitational waves can escape rapidly, leading to linearly 
increasing |a;i|. On the other had, as shown in the present paper, it is the "singularities" in the density profile of NSs 
that give rise to the above-mentioned distinguishing features in their w-mode spectra. 
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Appendix A: High-order HFA 



From p.8p we can expand F in power series of l/cj, i.e., 

p ^ ^(1) +^(2) _^ ^(3) ^ 



(Al) 
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where F<^''^ - 0{S/uj'') with k = 1,2,3,.... For purpose of reference, the lowest four orders of F'^'^^ are hsted as 
follows: 



^g-(.+A)/2 

Up- dr* ' 
m + Anr^p + {n + l)r 



g(t/-A)/2 r 4^ 



— + 87rrp - 



e^{n + 1) + e 



xm — Anr p 



dS_ 



— (n + 1) + ^ ^ 



2 + e-^(n + 1) (47rH(p + 2p) - m) 



2m 

+3(n+ 1) 



S. 



(A2) 
(A3) 
(A4) 

(A5) 



Substituting the expansion of F back into p.4p . we find a single second-order ODE for polar w-mode oscillation 
inside the star: 



drl 



+ (w^ -v^)s- + + + ...)= , 



(A6) 



where 



Vsf = 



4e^ 



(to + inpr-^)^ 
r — 2m 



47r/9r — 3m 



(A7) 



We have carried out numerical study to check the validity of (jA6l) and found that the second-order approximation 
(i.e., F = F^-^^f + F'^^'>) can readily reproduce polar w-mode QNMs with good accuracy and works for both leading 
and non-leading modes. As shown in Table XI, the accuracy of the perturbation scheme improves with increasing 
frequency and increasing order of expansion as well. 



Appendix B: Location of the maximum of dp/dr* 

Following directly from the TOV equations [H, [sl] , it can be shown that 



dp Ap + pf 



where 



h(r) = 1 dr m + Anr^p 

r(r) dr^ r{r - 2m) ^ ' 

is a smooth function in comparison with the term {p + p)'^ /p near the surface of a NS. Ignoring terms proportional to 
dh/dr^,, we find that 



dri p dr^ 



h{ry—^^[2-T{r)]. (B3) 



Therefore, the maximum value of |fip/fir*| is attained around the point r, — R2 where F = 2. 
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TABLE I: The first ten polar w-mode QNMs of a polytropic star witfi EOS p = lOOp'^ and a compactness 0.297. The first two 
columns are the exact values of lui and oji. The third and fourth columns are the approximate values of uuv and lui obtained 
from HEA. All uj are measured in units of . The last two columns are the relative percentage errors in uj^ and uii. 



exact uji- 


exact \(jJi\ 


TTT7A / 1 


-tlrA |Wi| 


lA, , // . 1 

\ouji-/uji-\ 


\ouji/uji\ 


0.5594 


0.3843 


0.6346 


0.4642 


13.45 


20.79 


0.4707 


0.0562 


0.4709 


0.0282 


0.04 


49.83 


0.6534 


0.1638 


0.6385 


0.1248 


2.29 


23.82 


0.8910 


0.2269 


0.8663 


0.2084 


2.77 


8.13 


1.1270 


0.2614 


1.1081 


0.2527 


1.67 


3.33 


1.3624 


0.2866 


1.3475 


0.2814 


1.10 


1.82 


1.5982 


0.3068 


1.5858 


0.3033 


0.78 


1.13 


1.8343 


0.3237 


1.8237 


0.3212 


0.58 


0.77 


2.0706 


0.3383 


2.0613 


0.3364 


0.45 


0.56 


2.3071 


0.3511 


2.2988 


0.3497 


0.36 


0.41 



TABLE II: The first ten polar w-mode QNMs of a realistic star with APR2 EOS and a compactness 0.196. The numerical 
results of u are given in units of M~^. The first two columns are the exacts values of uji and uji. The third and fourth columns 
are the approximate values of uji and oji obtained from HFA. All cu axe measured in units of . The last two columns are 
the relative percentage errors in Ui and uji. 



exact Wj- 


exact |wi| 


HFA CO, 


HFA Iwil 


|(5cJr/tJr| 


\SuJi/LL>i\ 


0.5050 


0.2971 


0.4575 


0.2064 


9.41 


30.53 


0.9209 


0.3826 


0.8827 


0.3676 


4.15 


3.92 


1.3067 


0.4378 


1.2824 


0.4319 


1.86 


1.35 


1.6840 


0.4803 


1.6660 


0.4773 


1.07 


0.62 


2.0565 


0.5139 


2.0421 


0.5123 


0.70 


0.31 


2.4266 


0.5405 


2.4146 


0.5396 


0.49 


0.17 


2.7966 


0.5613 


2.7862 


0.5608 


0.37 


0.09 


3.1681 


0.5786 


3.1589 


0.5782 


0.29 


0.07 


3.5413 


0.5945 


3.5330 


0.5942 


0.23 


0.05 


3.9155 


0.6106 


3.9081 


0.6102 


0.19 


0.07 
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TABLE III: Relevant physical parameters and values of b of model stars considered in the present paper. Here QS denotes the 
MIT bag model p ~ {p — 4_B)/3, where B = 56MeV fm""^ is used throughout the paper. C = M/R is the compactness of the 
star. The b values are calculated with I — 2. 





EOS Af 


C 


i?, R M 


n 


6'f, (Axial) 


6lf, (Polar) 




QS 


0.2644 24.617 11.292 2.986 4.3878 x lO"'' 


-2tt 


-37r 


p 


= lOp^ 1 


0.201 


5.354 2.873 0.577 


3.28 X 10~3 


-2.57r 


-3.57r 


p - 


= 5p3/2 2 


0.156 


33.602 16.961 2.646 


4.07 X 10~^ 


-37r 


-47r 


p = 


= 10p5/3 1.5 


0.195 


19.277 9.147 1.784 


2.24 X 10-5 


-2.757r 


-3.757r 


p = 


= lOp"/'^ 1.75 


0.195 


35.329 14.336 2.796 


1.50 X 10-^ 


-2.8757r 


-3.8757r 


p = 


= 5pi3/9 2.25 


0.144 


53.655 23.469 3.380 


2.38 X 10-^ 


-3.1257r 


-4.1257r 


p 


= p^^ 2.5 


0.117 


13.323 6.894 0.808 


2.32 X 10-*^ 


-3.257r 


-4.257r 



TABLE IV: Theoretical and numerical values of ci/tt, cq/tt in (|4.17p (l = 2) for the model stars tabulated in Table Hill The 
number enclosed in bracket is the value of co /n obtained directly from (|4.17|) using the theoretical value of ci and the frequency 
of the 20-th QNM. 





EOS 


method 


ci/7r(Axial) ci/7r(Polar) co/7r(Axial) co/7r(Polar) 




QS 


theoretical 


1 


1 


1 1.5 






numerical 


1.001 


1.001 


0.961(0.981) 1.462(1.481) 


p 


= lOp^ 


theoretical 


1 


1 


1.25 1.75 






numerical 


1.002 


1.002 


1.164(1.207) 1.665(1.708) 


p 


= 5p^/2 


theoretical 


1 


1 


1.5 2 






numerical 


1.002 


1.002 


1.410(1.448) 1.911(1.949) 


p = 


= lOp'/-' 


theoretical 


1 


1 


1.375 1.875 






numerical 


1.002 


1.002 


1.293(1.331) 1.795(1.831) 


p = 


lOp"/^ 


theoretical 


1 


1 


1.4375 1.9375 






numerical 


1.001 


1.001 


1.387(1.405) 1.887(1.905) 


p -- 


= 5pl3/9 


theoretical 


1 


1 


1.5625 2.0625 






numerical 


1.000 


1.000 


1.541(1.535) 2.040(2.035) 


p 




theoretical 


1 


1 


1.625 2.125 






numerical 


1.001 


1.001 


1.563(1.583) 2.063(2.083) 



TABLE V: Theoretical and numerical values of di, do in (|4.18|) (/ = 2) for the model stars tabulated in Table HITl The number 
enclosed in bracket is the value of do obtained directly from (|4.18[) using the theoretical value of di and the frequency of the 
20-th QNM. 





EOS 


method 


c?i (Axial) c?i (Polar) 


do (Axial) 


do (Polar) 




QS 


theoretical 


-1 


-1 


-0.662 


-0.662 






numerical 


-0.997 


-0.997 


-0.677(-0.664) 


-0.678(-0.664) 


p 


= lOp^ 


theoretical 


-1.5 


-1.5 


-0.343 


-0.343 






numerical 


-1.490 


-1.490 


-0.387(-0.347) 


-0.388(-0.347) 


p 


= 5p-V2 


theoretical 


-2 


-2 


-0.328 


-0.328 






numerical 


-2.003 


-2.002 


-0.320(-0.330) 


-0.320(-0.330) 


p -- 


= 10p5/3 


theoretical 


-1.75 


-1.75 


-0.175 


-0.175 






numerical 


-1.749 


-1.748 


-0.183(-0.177) 


-0.183(-0.177) 


p = 


lOp"/^ 


theoretical 


-1.875 


-1.875 


-0.020 


-0.020 






numerical 


-1.890 


-1.892 


0.045(-0.018) 


0.052(-0.018) 


p -- 


= 5p"/9 


theoretical 


-2.125 


-2.125 


-0.313 


-0.313 






numerical 


-2.164 


-2.166 


-0.138(-0.303) 


-0.129(-0.302) 


p 




theoretical 


-2.25 


-2.25 


-0.662 


-0.662 






numerical 


-2.281 


-2.279 


-0.525(-0.655) 


-0.533(-0.654) 
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TABLE VI: The axial w-mode QNMs {j = 3 to j = 12) of a reahstic NS (with EOS A compactness of 0.27 and central 
density pc = 2.227 x 10"^) are tabulated. The first two columns are the exacts values of ujr and uji. The third and fourth 
columns are the approximate values of cjr and Ui obtained from the the unified potential in (|5.2p . All uj are measured in units 
of Af~^. The last two columns are the relative percentage errors in uji and uJi, respectively. 



exact 


exact uj[ 




\uj\ 1 

1 1 


1 SlJ-,- I UJy 1 


\SuJi /iJi 1 


1.4031 


0.3690 


1.4105 


0.3694 


0.52 


0.13 


1.7543 


0.3974 


1.7604 


0.3973 


0.35 


0.03 


2.1065 


0.4205 


2.1117 


0.4202 


0.25 


0.07 


2.4593 


0.4407 


2.4639 


0.4404 


0.18 


0.07 


2.8127 


0.4593 


2.8167 


0.4589 


0.14 


0.07 


3.1662 


0.4772 


3.1698 


0.4769 


0.11 


0.07 


3.5195 


0.4949 


3.5228 


0.4946 


0.09 


0.06 


3.8724 


0.5128 


3.8754 


0.5125 


0.08 


0.06 


4.2245 


0.5313 


4.2273 


0.5310 


0.07 


0.06 


4.5752 


0.5501 


4.5778 


0.5499 


0.06 


0.05 



TABLE VII: The polar w-mode QNMs {j = 3 to j = 12) of a realistic NS (with EOS A compactness of 0.27 and central 
density pc — 2.227 x 10~^) are tabulated. The first two columns are the exacts values of uir and ua. The third and fourth 
columns are the approximate values of uui and cui obtained from the the unified potential in (|5.2p . All oj are measured in units 
of . The last two columns are the relative percentage errors in ujt and coi, respectively. 



exact ujr 


exact \uji\ 




\uJi\ 


\6uJr/uJr\ 


\6uji/uji\ 


1.2266 


0.3507 


1.2142 


0.3381 


1.01 


3.59 


1.5782 


0.3837 


1.5695 


0.3764 


0.55 


1.90 


1.9302 


0.4093 


1.9235 


0.4044 


0.35 


1.18 


2.2828 


0.4308 


2.2774 


0.4273 


0.24 


0.81 


2.6359 


0.4500 


2.6314 


0.4474 


0.17 


0.59 


2.9894 


0.4682 


2.9856 


0.4661 


0.13 


0.45 


3.3429 


0.4860 


3.3395 


0.4843 


0.10 


0.35 


3.6960 


0.5038 


3.6930 


0.5024 


0.08 


0.29 


4.0486 


0.5220 


4.0459 


0.5208 


0.07 


0.24 


4.4001 


0.5407 


4.3977 


0.5396 


0.05 


0.20 
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TABLE VIII: The values of Rm/R*, R2/R* and R^/R^ for the two realistic stars constructed with the analytic approximation 
of the SLy EOS [5^ with compactness C — 0.306, 0.228 are tabulated. The value of R, is obtained from from the best linear 
fit to the data shown in Fig. 2(a) 



c 


Rrn 1 R* R2 1 R* R* 1 R* 


0.306 
0.228 


0.980 0.978 0.984 
0.957 0.950 0.958 



TABLE IX: Evaluation of Z = 2 polar w-modes from the corresponding axial w-modes for a model C NS (with compactness 
C = 0.186 and central density pc = 8.752 x 10~*). The first two columns show the exact values of ujr and uji. The third and the 
fourth columns are respectively approximate values of ojr and Ui obtained from (|2.8p . The last two columns are the approximate 
values of ojr and ua obtained from corresponding axial w-modes with I — 3. All lj are measured in units of A4~^. 



exact UJr 


exact 








UJr 






0.9135 


0.2286 


0.9134 





2283 


0.9118 





2287 


1.0796 


0.2398 


1.0795 





2396 


1.0781 





2399 


1.2453 


0.2493 


1.2453 





2491 


1.2440 





2494 


1.4110 


0.2574 


1.4110 





2573 


1.4099 





2575 


1.5769 


0.2645 


1.5769 





2644 


1.5758 





2646 


1.7429 


0.2711 


1.7430 





2710 


1.7420 





2711 


1.9091 


0.2773 


1.9091 





2772 


1.9083 





2773 


2.0754 


0.2831 


2.0754 





2831 


2.0746 





2832 


2.2419 


0.2888 


2.2419 





2887 


2.2411 





2888 


2.4085 


0.2943 


2.4086 





2943 


2.4078 





2943 
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3 3.5 4 4.5 

(b) 

FIG. 1: The axial (circles) and the polar (crosses) w-mode QNMs of the polytropic stars considered in Table [Till are shown, 
(a) _R*cjr — an/2 is plotted against the mode number j, j = 1, 2, 3, 20. For clarity the data of each polytropic star are 
upward-shifted by 5, 10, 15, 20, 25, 30 for Af = 1,1.5,1.75,2.0,2.25,2.5, respectively, (b) is plotted against \n{R,ujr) 

for 5 < j < 20. The solid straight lines in (a) and (b) are the corresponding theoretical lines predicted by (|4.14p and 1)4. 15|) . 
respectively. 
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(b) 

FIG. 2: The axial (circles) and the polar (crosses) w-mode QNMs of two realistic stars obeying the analytic approximation of 
the SLy EOS proposed in Ref. 59] with compactness C — 0.306, 0.228 are shown, (a) i?,tJr — i7n/2 is plotted against the mode 
number j, j — 1, 2, 20. For clarity for the case C = 0.306 the data points are all shifted upward by 5. The two straight solid 
lines shown in the figure are the best linear fits to the two sets of data, respectively, (b) |-R*a;i| is plotted against ln(J?»ajr) for 
5 < j < 20. 
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FIG. 3: dp/dr, (normalized by pcjR*, with pc being the density at the center of the star) is plotted against r* (normalized by 
i?*) for two realistic stars obeying the analytic SLy EOS proposed in Ref. [s^ with compactness C — 0.306, 0.228. As labelled 
in the graphs, \dpldr,\ attains maximum at r* = 0.98i?, (0.957i?,) for the case C = 0.306 (0.228). 




0.4 

3 3.2 3.4 3.6 3.8 4 4.2 
ln(i?*6Jr) 



FIG. 4: A/", Mm, and TVe are plotted against ln(_Rtajr) for two SLy stars with C 



= 0.228, 0.306, respectively. 



